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Nanoscale hydrogen particles in superfluid helium track the motions of quantized vor- 
tices. This provides a way to visualize turbulence in the superfluid. Here, we trace the 
evolution of the hydrogen from a gas to frozen particles migrating toward the cores of 
quantized vortices. Not only are the intervening processes interesting in their own right, 
but understanding them better leads to more revealing experiments. 
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I. INTRODUCTION 

What makes liquid helium particularly useful for stud- 
ies of fluid dynamics and turbulence is its sma ll viscosity 
(jNiemela and Sreenivasanl . 120061: ISkrbekl . 120041 ) . Another 
interesting aspect of liquid helium is that in the super- 
fluid phase circulation is quantized, and turbulence arises 
in it thro ugh the formatio n and interaction of quantized 
vortices ( Fevnmanl . 1955 ). The superfluid is present in 
liquid helium-4 at temperatures between OK and T\, 
w here T\ 2.1 7K when the liquid and vapor coexist 



(|Donnellvi . 119911) . 

The strides made in understanding classical turbulence 
have not been matched in quantum turbulence. Chal- 
lenges lie in both the theoretical description of the fluid 
motions and in t he acquisition and interpr etation of ex- 
perimental data (IVinen and Niemelal . [2002h . In an effort 
to surmount the latter hurdle, one of the authors (GB) 
developed a novel hydrogen-particle technique for visual- 
izing the internal motions of superfluid helium. The tech- 
nique involves the generation of the pa rticles in situ, and 
sequential recording of their positions ( Bewlev . 20091) . 

What experimenters observe is the following. 
Nanoscale hydrogen particles appear in liquid he- 
lium close to its boiling point, after injection of a 
mixture of hydrogen and helium gas through a small 
tube. When the temperature is lower than T\, the 
particles gravitate toward quantized vortex cores. The 
entrapment does not occur when the temperature is 
close to T >, bu t only at significantly lower temperatures 



(jBewlevi . 120061) . 

The motions of quantized vortices become visi- 
ble because particle s beco me trapped on their cores 
([Parks and DonnelM I1966I ). Once they are visible, it is 



possible to observe t he vortex lattice in r otating helium 
( Bewlev et al. , 20061) , vortex re connection ( Bewlev et al. , 
2008 ). and vortex ring decay (jBewlev and Sreenivasan . 
2009J). An increasing number of laboratories around the 



world are employing similar techniques. They do so for 
techno logical reasons , such as for the generation of nano- 



wires ([Gordon et all 120121 ) , but above all to study the 



fluid dynamics of helium. 

Since the technique is employed in many ways, 
it makes sense to better understand it. The self- 
organization of hydrogen onto quantized vortices involves 
non-equilibrium processes occurring on many length 
scales. These start from the microscopic physics of 
demixing and freezing that produce the particles, to the 
mixing of the particles into the liquid and the migration 
of the particles toward the cores of vortices, to the macro- 
scopic turbulent motions of the helium itself. This paper 
proposes descriptions for each process. These should be 
understood to be conjectural, as we have not tested the 
ideas experimentally. 

Specifically, we set out to explain how the hydrogen 
particles form, how they end up in the liquid rather than 
the gas phase of helium, and why vortices do not entrap 
them when the temperature is close to T^. In Sec. III. Al 
we describe the experimental apparatus. In Sec. III.B1 we 
survey briefly recent work that exploits particle-imaging 
to study superfluid helium dynamics. In Sec. IIII.AI we 
describe the phase transitions that produce hydrogen 
particles, and in Sec. IIII.BI the process dispersing the 
particles in the liquid helium. In Sec. IIII.C1 we argue 
that thermodynamics does not limit the trapping of par- 
ticles unless the temperature is very close to T\, and in 
Sec. IIII.DI we argue that particle dynamics, rather than 
thermodynamics, limit the ability to visualize vortices. 



Finally, in Sec. II VI we comment on these interpretations. 



II. HYDROGEN PARTICLE GENERATION REVISITED 



A. Experimental setup 



La Mantia et ali (j2012ri and iBewlevI (|2009h describe 



the relevant experimental apparatus, which we briefly re- 
view here. Hydrogen particles are created from a room 
temperature mixture of hydrogen gas diluted with he- 
lium. The method works also with other gases, such as 
deuterium, in place of hydrogen. The hydrogen- helium 
mixture ratio controls the size of the particles, the size 
increasing with the ratio of hydrogen to helium. Mixture 
ratios of the order of 1 H2 : 100 He by volume produce 
particles about 1 /im in diameter. The mixture is re- 
leased into liquid helium through valves that control the 
injection rate and duration. For a 600 mL volume of liq- 
uid helium in a cryostat, 20 mg of the mixture injected 
over a period of 5 s yields a final particle volume fraction 
about 10 5 times smaller than the total volume, which 
corresponds to a few million particles per cm 3 . The gas 
mixture passes through a tube that opens somewhere be- 
neath the free surface of the helium, so that the mixture 
bubbles through the liquid. 



process (jBewlev et all 120081 ) , while at high number den- 
sities, networks are for med where cr o ssing vortices stick 
rather than reconnect ( Bewlev et ali . [20061 ) . 

The behavior of particles in supcrfluid helium and in 
particular the interaction between particles and quan- 
tized vortices have been the subject of many numerical 
and theoret ica l inve stigations. Notable examples include 
Poole et alx s ([2005) discussion of the motions of parti- 
cles in the context of the equations of particle motion. 
Fujiyama et ali (|2007l ) demonstrates the pinning of quan- 
tized vort ices by solid spheres, a n interaction explored 
furthe r bv lKivotides et al. (120081). Barenghi and Sergeev 
(2009) and iBewlev and Sreenivasa nl (<2009h consider the 
influence of bound particles on quantized vortex dynam- 
ics. This work inspired interpretations of experimental 
observations. It also argued that particles are trapped 
when the relative velocity between the normal fluid and 
the quantized vortex is not too large, and that the pro- 
cess of particle trapping by a quantized vortex curves 
the vortex core and causes Kelvin waves along the core. 
Once trapped, the particles increase the strength of the 
coupling between the quantized vortices and the normal 
fluid. 



III. PHYSICAL MECHANISMS 



B. Recent history 

We cite a few p apers to give an im pression of 
the field's breadth. iPaoletti et ali ([20080 shows that 



the distribution of particle velocities in superfluid he- 
lium is non-Gaussian, in contrast to the near Gaus- 
sian statistics of particle moti ons in classical flows. 
Chagovets and Van Sciverl (|201ll ) measures in a counter- 
flow the heat flux required to dis l odge hydrogen parti- 
cles from vortices. Jin and Maris ( 201 lh describes some 
interestin g behavior of dus t part icles in supcrfluid he- 
lium, and La Mantia et al. ( 20121 ) confirms that hydro- 
gen particle velo cities equal the norm al fluid velocity at 
low heat fluxes. Fuzier et al. (120081) shows that super- 
fluid helium can be used to segregate particles by size. 
Quo et ali d200flh images electron bubbles and speculates 
that this technique could be used to track quantized vor- 
tices. Metastable helium molecules c a n also be tracked, 
an approach employed by Guo et al. ( 2010h to measure 
the velocity profile of superfluid helium in a channel. 
Gordon et all (|2012l) discusses the role of vortices in the 
generation of metallic nanocrystals and nanowircs in su- 
pcrfluid helium. Taken together, the collected works in- 
dicate that under some conditions, such as in weak coun- 
terflows, quantized vortices trap particles on their cores, 
while in other conditions, the particle- vortex interactions 
yield more complicated behavior. At low particle number 
densities and when the particles are trapped, the particle 
motions reveal that vortex rcconncction is a self-similar 



No one has yet found a way to observe particle forma- 
tion directly. In the following, we lay out our personal 
view of the underlying mechanisms, which involve phys- 
ical problems worth studying in their own right. 



A. Particle generation 

In a typical hydrogen particle generation experiment, 
gas is injected at a rate of J = 4 x 10~ 3 g/s through 
a pipe with a diameter of about 0.15 cm. Given the 
mass density, p g w 1.5 x 10 _3 g/cm 3 , this amounts to 
a flow velocity of 14cm/s. For aim tall cryostat, the 
gas takes about 7 s to exit the tube into the liquid he- 
lium bath. As it descends through the injector into 
the cryostat, the mixture of hydrogen and helium gas 
has a higher temperature than the ambient helium bath. 
The generation of particles should therefore be under- 
stood as a problem of cooling predominantly helium gas 
bubbles that cool until hydrogen condenses. In such 
a setting, we expect that the initial h ydrogen particle s 
are created by homogeneous nucleation (|Doremusl . 119851) . 
This process occurs in analogous binary fluids cooled 
into t he miscibility gap ( Lapp et al . 20121 Vollmer et al. . 
20071 ). At the low temperatures of helium, quantum 



tunneling may promote cry stallization of the hydrogen 
( Levi and Mazzarerloll200lj ). 

Specifically, the hydrogen starts as a warm gas strongly 
diluted by helium. After injection, the gas mixture bub- 
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FIG. 1 Using a schematic phase diagram, we follow the path 
of a helium bubble as it cools while rising through liquid he- 
lium. When the influence of the hydrogen is neglected, the 
liquid-vapor coexistence curve marks the boundary between 
the gaseous (below) and liquid (above) phases of helium. The 
A-line marks the boundary between normal liquid helium (to 
the right), and the liquid phase with superfluid properties (to 
the left). The pressure difference between the bath and the 
bubble, Ap = 2a /R 2 x 1 0~ 6 bar, is due to surface tension 
(|Atkins and de Paulal . [20021 . sec. 6.9), and has some time de- 
pendence that we do not show. Here, a is the surface tension 
and R is the radius of the bubble. When the bubble liquifies, 
hydrogen particles suspended in the bubble are released. 



bles cool as they rise through the liquid helium. The 
melting and boiling temperatures of hydrogen at the rel- 
evant pressures (e.g. 14 K and 22 K at lbar) exceed the 
boiling point of helium (4.2 K at lbar), so that the hy- 
drogen condenses out of the mixture even while the he- 
lium fraction remains gaseous. Cooling proceeds even 
after particle nucleation. In this way, the condensation 
is similar to the formation and growth of rain droplets 
in orographic rain formation, w here c l ouds cools slowly 
due to adiabatic expansion (e.g. lShawl . [200l . to the for- 
mation of droplets in laboratory experiments where bi- 
nary fluid mixtures are slowly pushed into a phase coexis- 



tence region by an applied temperature drift ([Lapp et al 
2012HVollmer et al\ . \200lb . and to the formation of slow 



dust in Enceladus' plume from the condensation of wa- 



ter ( Schmidt et al 



2008). Understanding of phase sep- 



aration in such settings is presently a subject of active 
research due to the importance of clouds in climate dy - 
namics ( Bodenschatz et al. , 2010l ; IPierrehumbertl 2010() . 
Based on the studies of phase separation in the pres- 
ence of a temperature dr ift, we expect a w ide distribu- 



tion of particle sizes (e.g. lLapp et all 120121) . The distri 



bution results from the interplay o f growth and ripenin g 



(jLifshitz and Slvozovl[l96l1:ISlezovl . [20091: IWagnei" 



B. Deposition of the particles into the liquid 

From the analogy with boiling water, one might ex- 
pect that the bubbles would grow in size and then break 
through the liquid helium surface, dispersing the hydro- 
gen particles in the gaseous phase of helium. However, a 
bubble is not stable in a liquid with a free surface, due 
to the high free energy associated with the surface ten- 
sion. In experiments, bubbles grow smaller at a rate that 
is limited by how quickly the heat of the bubble can be 
transferred from it to the surrounding liquid. This oc- 
curs even if the system is on the liquid- vapor coexistence 
curve. 

The bubbles cool at approximately constant ambient 
pressure, as illustrated in Fig.[TJ due to the low density of 
liquid helium. One must consider six contributions to the 
heat of the bubble. These are due to the temperature dif- 
ferences between the surrounding helium and the hydro- 
gen and helium in the bubble, the heat of condensation 
and of freezing of the hydrogen, the heat of condensation 
of the helium, and the heat released as the surface area 
of the bubble decreases. 

The key point is that each contribution to the free 
energy of the bubble constitutes an amount of heat 
that must be transfer r ed to the surrounding liquid 
( Chaikin and Lubenskv . 2000h . The rate at which the 



bubble collapses is determined by the rate at which this 
heat can be carried away from the bubble by diffusion or 
convection. Initially, the mass of the bubble will increase 
due to this heat transfer, because some of the surround- 
ing helium will evaporate into the bubble. Eventually, 
however, the temperature of the bubble is low enough 
that it condenses. Here, we consider each contribution in 
turn. 

The heat due to the initial warmth of the droplets is 

AQ heat ps — R 3 Pg C p AT > 4 x lO" 4 J 

where p g w 1.4 kg/m 3 , C p « 1.2 x 10 2 m 2 /s 2 K, and 
AT > 10 K, are the gas density, its heat capacity and 
a conservative estimate of the initial temperature differ- 
ence, respectively. 

The latent heat of the phase transitions include those 
from the fluidization and solidification of hydrogen, and 
from the fluidization of the helium in the bubble. They 
A^ idize = 4.5 x 10 2 J/kg, Aq^ idiiy = 6-0 x 10 4 J/kg, 
and Aq^ idizc = 2.1 x 10 4 J/kg, respectively. For the 
cooling of the bubbles they give rise to the contributions 
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where p, <fi and Aq denote the material constants for the 
respective materials and transitions; in particular, <p 
1 200 is the mass fraction of hydrogen in the droplets. 
It is so small that cf> 1 for the helium. The latent 
heat released upon going through these phase changes, 

AQcondcnsation = ^Qfluidizo + ^^solidify + ^QfhfidizC 

accounts for about half of the total heat, AQ t ot = 
AQcondensation+AQheat, to be released from the droplets. 

The bubbles collapse when, upon rising, they dissipate 
the excess heat AQtot into the fluid helium bath before 
they reach the interface to the macroscopic gas phase. 
Note that the heat from the temperature mismatch of 
the bubbles exceeds the contribution from the surface 
energy by five orders of magnitude, such that the latter 
may be neglected. 

We estimate how long it takes for the bubble to col- 
la pse by comp a ring to the experimental observations 



of Maris et al. (1 1 98 7 ) and to the numerical study of 



Unlusu et al. ( 20081 ). Both studie s treated pure hydro gen 



droplets in a liquid helium bath. iMaris etali ( 19871 see 



their Fig. 3) observed that the te mperature of d r oplets 



equilibriated in about 0.2 seconds. lUnlusu et all (|2008f ) 
found that the droplets released their heat in a few sec- 
onds. We must bear in mind that the heat content of the 
hydrogen droplets is more than 25 times that of helium 
bubbles, due to the very high heat of solidification of the 
hydrogen. Our mixtures of He and H2 will therefore cool 
faster. 

How long does it take for the bubbles to reach the free 
surface in the cryostat, setting aside any changes in size? 
Due to the low viscosity of liquid helium, the bubbles 
generate a turbulent wake upon rising. The resulting 
rising velocity is approximately 
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One can check that the assumption of a turbulent wake 
is justified. Given the density, p/ = 140 kg/m 3 , and the 
viscosity, 77 = 4 x 10 -6 Kg/ms, of liquid helium, the ris- 
ing bubble has a Reynolds number Re = pfUR/rj « 500, 
which lies well in the turbulent regime. Consequently, 
when they are released at a depth of about ten centime- 
ters, the bubbles need sa 10 seconds to rise to the surface. 
This is considerably longer than the time needed for the 
heat content of the bubbles to be carried away. 

We conclude that the helium bubbles condense and re- 
lease their hydrogen particles to the liquid helium before 
they re ach the t op of the liquid helium column. Accord- 
ing to (jSilveral Il984h . the hydrogen particles have no 
difficulty crossing the liquid-gas interface of helium, so 
that they will move to the liquid phase continuously as 
the bubbles shrink, rather than being caked into a sin- 
gle large aggregate before release. Hence, we observe a 
much larger number of small hydrogen clusters than the 
number of bubbles released into the cryostat. 
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FIG. 2 As a consequence of the rapid rise of the condensation 
energy e, condensation arises in a very narrow temperature 
interval. According to our model, Eq. (|A4|I . which addresses 
the thermodynamic equilibrium distribution, the vortices will 
be densely packed less than 50/iK below T\. 



C. Thermodynamics of particle trapping 

We have already observed experimentally that quan- 
tized vortices trap hydrogen particles only when the tem- 
perature of the helium is sufficiently far below T\. We 
now show that thermodynamics does not explain this 
link. The strength of the interaction between quantized 
vortices and particles is proportional to the mass density 
of the superfluid part of the helium, which rises sharply 
from zero as temperatures fall below the phase transition 
temperature. However, for some range of temperatures 
immediately below T\, the interaction is still negligible 
relative to the thermal energy, fcgT, so that particles re- 
side on the vortex cores only accidentally. It follows that 
as the temperature falls, the transition between free par- 
ticles and trapped particles occurs later than the transi- 
tion from normal to superfluid helium. 

In the absence of other effects, such as macroscopic 
flows, the nature of the trapping transition can be de- 
termined with a lattice model where the lattice con- 
stant amounts to the characteristic size of the hydro- 
gen particles. In this context, we say that vortices run 
through a fraction <j>y of the lattice sites, such that 
there are Ny = <j>yN sites occupied by vortices, and 
Nb = (1 — 4>v)N sites making up the fluid bulk, where 
N is the total number of sites in the lattice. In the Ap- 
pendix, we show that the fraction of particles on vortices 
is then 



R 



4>v 



4> H + exp(/3e) ' 



(1) 



where <pu is the volume fraction of particles, f3 = 1//cbT, 
and e is the change in the system energy due to the trap- 
ping of a single particle on a vortex (which is negative). 
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The factor e contains the most important temperature 
dependence, which is that due to the ra pid change in den- 
sity o f the superfluid part near T\ fe.g. lTam and Ahlers . 
19871 ). It can be seen in Fig. [2] that according to this 
model, the vortices saturate with particles at tempera- 
tures only 50 pK below T^. Because experimentally we 
observe little particle entrapment up to 20 mK below T A , 
it follows that the ability to track vortices is not limited 
by thermodynamic considerations. Instead, it is domi- 
nated by the dynamics of particle trapping. 



D. Dynamics of particle trapping 

At temperatures far enough below the phase-transition 
temperature, q uantized vortices begin to gather particles 



on th eir cores (jBewlev et all l2006t iParks and Donnelly , 
1966(), provided the relative velocity between the quan- 
tized vortices and the normal fluid is not too high 



(IChagovets and Van Sciverl . 120111 ) . We estimate here the 
time required for a newly-formed length of quantized vor- 
tex to collect enough particles to make it visible. 

Consider a stationary vortex that extends through 
a fluid filled with initially stationary particles, whose 
total volume fraction relative to the fluid is tpn, and 
whose characteristic diamter is d. The distance between 

1 /3 

particles in the fluid is then typically l a = d/{2(j>^ ). 
At a large distance r from the core, the force exerted 
by the vortex on the particle is approximately F v = 
— (l/6)nd 3 d r p, where p = — p s n 2 /8ir 2 r 2 is the pres- 
sure, p s is the density of the superfluid, n = h/m ~ 
1 x 10~ 3 cm 2 /s is the circulation of the vortex, h is 
Planck's constant, and m is the mass of a helium-4 atom. 
Let the vortex be considered marked when all parti- 
cles within this distance from the vortex centcrline have 
reached the core. It can be readily confirmed that this 
corresponds to a core with particles every 20 diameters, 
for a typical particle volume fraction of 10~ 5 . If the par- 
ticle's inertia is small, we can balance the vortex pressure 
gradient force by the drag force along its whole trajec- 
tory. In this picture, the particles are drawn toward the 
core by the superfluid pressure field, but are damped by 
their motion relative to the stationary, viscous normal 
fluid. The drag force for small particles was given by 
Stokes as Fs = Zitpdv, where p is the viscosity of the 
normal fluid and v is the veloci t y of t he particle relative 
to the normal fluid (|Batchelorl 119671 ). The velocity, v, 
is induced by the pressure gradient such that F v = F$. 
Altogether, the time scale for marking is given by inte- 
grating the velocity from a distance l to the vortex core, 



(2) 



The time scale is temperature dependent through the ma- 
terial properties, and particularly through p s , which rises 





FIG. 3 The figure shows the dependence of the vortex mark- 
ing time, Tm, given by equation [2] on the temperature. Here, 
the particle diameter, d, is 0.5, 1, and 2 microns for the lower, 
middle and upper curves, respectively, and the particle vol- 
ume fraction, (j>H, is lxlO -5 . For lower volume fractions, the 
marking time is longer. The inset shows the same curves on a 
log-log scale, and as a function of the temperature difference, 
T A - T. 



from zero for temperatures decreasing from the phase 
transition temperature, T A . Figure [3] shows the depen- 
dence of r m on temperature for values of d and <pjj typi- 
cal in an experiment. The sharp decrease in the trapping 
time as the temperature decreases may explain the ob- 
servation that marked vortices typically do not appear 
within about 20 mK of the phase transition temperature. 
That is, they might become visible if the experimenter 
waited long enough. 



IV. DISCUSSION AND CONCLUSION 

The journey of hydrogen from a molecular gas to crys- 
tallites that are trapped on quantized vortices encom- 
passes a hierarchy of problems on many length scales. 
Novel problems arise in applied mathematics through 
the analysis of nonlinear evolution equations, in non- 
cqulibrium statistical physics through the condensation 
of hydrogen and the collapse of bubbles, and in fluid 
dynamics through the mixing of particles into a fluid, 
the trapping of particles onto vortices, and the turbu- 
lent motions of the vortices. These challenges call for 
a combined effort from experimentalists and computa- 
tional physicists to provide data, and from theoreticians 
to identify approaches and questions. 

In this spirit of providing guidance for future exper- 
iments, where efficient particle generation is desired or 
where manipulation of the particle sizes and densities 
are necessary, we make the following propositions. There 
should be many small injectors, which produce small 
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bubbles that disappear quickly. The injectors should be 
as far below the surface of the liquid helium as is pos- 
sible, and the injector tube could be coiled within the 
cryostat in order pre-cool the mixture it flows into the 
helium liquid. The cryostat should be as tall as possible. 
The particles should be as small as possible in order to 
become trapped quickly. 
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are 

lattice sites. Vortices run through a fraction <j> v of these 
sites such that there are N v = 4> V N sites occupied by a 
vortex and Nb = (1 — 4>v)N sites making up the fluid 
bulk. 

We choose a description in terms of a grand-canonical 
potential f2(/x,T) where T is the temperature and the 
chemical potential \i will be adjusted to achieve that a 
fraction <pH of the sites are occupied by particles. Finally, 
particles that condense onto a vortex change the system 
energy by an amount e < 0. 



Appendix A: Condensation of Hydrogen crystallites on 
vortices 

To support our discussion of the trapping of the crys- 
tallites on the vortices in Sec. IIII.C1 we estimate their 
equilibrium distribution based on a simple lattice model, 
where the trapping amounts to condensation into bound 
states of particles on the vortex. The binding energy is e, 
and the particles have a monodisperse size distribution 
with diameter d. We assume that the particles do not 
interact (in particular, they do not merge). 



1. Model 



2. Partition functions 

For non-interacting particles the grand-canonical po- 
tential f2(^t,T) is the sum of the functions Qb and fly 
describing particles in the bulk and on the vortices, re- 
spectively, 



n(ii,T) = n B (j i ,T) + n v ( Ji ,T), 



We consider a lattice model with a lattice constant of 
the order of the particle diameter d. In a volume V there 



They take the respective forms ( C alien . 1985 : Widoml . 
20021) 



Q B (/J>,T) = -k B T In 
fV(At,T) = —kBT In 



'N B 

E 

.71=0 

'N v 

E 

.n=0 



V 



= -fc B T(l-0„)— ln[l + e-^] 



V \ e -/3(H-e)« 



V 

-k B T(f) v — In 
d 6 



1 + e 



where ks is the Boltzmann constant. where 

The total number of particles, Nh, is obtained as the 
negative derivative of f2(//,T) with respect to fj,, such that 
the volume fraction of the particles amounts to 



M = cT^ L (A2b) 

N H d 3 = d anp,T) E = e -pe (A2c) 

V V d\x 

<j> v EM , {1-4> V )M 
Y^EM + l + M ( ' 
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4. Fraction of condensed particles 



400 




From Eq. (|A2a[) we see that the volume fraction of 
particles condensed onto the vortices amounts to 



V dp 



<t> v EM 
1 + EAI 



To find the fraction of particles condensed to the vortices 
we divide by the overall volume fraction of particles, 4>h- 
Together with Eq. (|A"3|) . 
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FIG. 4 The condensation energy — t in units of fcsT as a 
function of the temperature (c/. Eq. (|A5[) '). 



3. Fixing the chemical potential 

To compare to the experiments we have to fix the 
chemical potential, i.e., the parameter M in Eq. (|A2a[) . 
in order to fix the particle volume fraction to the partic- 
ular value 4>h- The exact expression for M(4>h) can be 



b H + exp(/3e) 

f <f>v j f° r "C 1 (accidental) 

[ (j)v/(f>H , for /3e < -1 (full) 



(A4a) 



(A4b) 



determined straight-forwardly by multiplying Eq. (IA2a[) 
by (1 + ME) (I + M) and solving the resulting quadratic 
equation. 

For the following discussion, however, we observe that 
<S <f)y <C 4>h <C 1 in order to find a more transparent 

(IA2al) in the form 



expression. Rewriting Eq. 



M 



>n = 



l + M 



rM 



E 



1 



1 + EM l + M 



This result states that for \f3e\ <C 1 a fraction cf>y of 
the particles reside on sites with vortices. In this limit 
the energy change e for residing on a vortex is small as 
compared to fcsT so that particles do not feel the vor- 
tices. They are randomly distributed, and the particles 
are found on a vortex only by coincidence. 

In the opposite limit j3e <C — 1 the particles are strongly 
attracted to the vortices. However, since <f>jj 3> 4>Vi not 
all particles can be placed on a vortex. Eventually, when 
all sites Ny on the vortices are occupied, the fraction 
of condensed particles amounts to the number of lattice 
sites on vortices, Ny = N(f>y, over the number of parti- 
cles, Nh = N(f>H, such that Ny/Nn = <t>v/<i>H- 



we see that the expression in square brackets is of order 5 ' Com P ariso " to the experiments 

(3e/(l + M) 2 for small \/3e\ where E ps 1. For cV < 1 one 



must then have a value of M of the order of and the 
second term in the sum may be neglected. Consequently, 



M 



The cross-over from freely-moving particles to con- 
densing particles arises when <f>H ~ e ,3e , i.e., for e sa 
ksT In <j)H- In order to pin down the locus of the transi- 
tion we observe that 



1 



te = 6x lir 6 



On the other hand, when (3e becomes much smaller than 
— 1 such that E ^> 1, the second term in the square 
bracket may be neglected and the one in the denominator 
of its first term will approach EM, such that 



L,y 2 

a 

V 



1.6 x 10 4 



io- 



1.6 x 10" 



M 



<i>H - <\>V 

1 - (f>H — <f>V 



Finally, the binding energy -e/^T is estimated as the 
amount of kinetic energy of superfluid which would still 
be circli ng around the vortex core in the absence of the 
droplet ( Parks and Donnelly . 19661) . 



The exact solution varies monotonously between these 
limits, and in view of <C e>y <C 4>h < 1 we may safely 
assume in the following that 



-iTdp s 



( h 



1 - 



a 

1 + ^ 



1/2 



asinh ( — 
a 



(A5) 



M 



(A3) 



where a k 10 cm is the diameter of the vortex core, 
"^Hc ~ 6.688 x 10 _24 g is the mass of helium atoms, and p s 



8 



is the mass density of the supcrfiuid com ponent of helium 
which is approximated according to Tarn and Ahlersl 
(|l987l ). Due to the rapid rise of the mass density of 
the supcrfiuid phase upon crossing T\ this function rises 
rapidly when T drops below Ta (Fig. 2]). Inserting this 
data into Eq. (|A4ap provides the prediction for R shown 
in Fig. dJ 
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